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Abstract:  Different  release  models  can  yield  significantly  different  dy¬ 
namic  concentration  profiles  in  a  room  depending  on  the  release  rates 
chosen.  Results  from  two  different  cases  demonstrated  significantly  differ¬ 
ent  concentration  profiles  in  the  room  of  interest.  This  work  was  under¬ 
taken  to:  (l)  Critically  review  the  current  model,  (2)  formulate  a  new,  1st 
order  Algebraic  model,  and  (3)  use  experimental  data  to  validate  the  mod¬ 
eled  theory.  This  report  documents  preliminary  work  that  was  suspended 
after  Fiscal  year  2005  (FY05).  No  follow-on  research  is  currently  funded. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes.  Citation 
of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products.  All  product 
names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 

DESTROY  THIS  REPORT  WHEN  NO  LONGER  NEEDED.  DO  NOT  RETURN  IT  TO  THE  ORIGINATOR. 
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Preface 

This  study  was  conducted  for  the  Engineer  Research  and  Development 
Center,  Construction  Engineering  Research  Laboratory  (ERDC/CERL), 
U.S.  Army  of  Corps  of  Engineers  under  6.1  Basic  Research,  “Resuspension 
Physics  of  Fine  Particles”  The  technical  monitor  was  Dr.  John  M.  Cullin- 
ane,  CEERD-EM-J.  This  study  was  also  supported  in  part  by  the  Project 
GDiKGD  “Bldg  Chem-Bio  Protection  Modeling  &  Simulation.”  The  techni¬ 
cal  monitor  was  Dr.  Paul  Howdyshell,  CEERD-CV-ZT. 

The  work  was  performed  by  the  Energy  Branch  (CF-E)  of  the  Facilities  Di¬ 
vision  (CF),  Construction  Engineering  Research  Laboratory  (CERL).  The 
CERL  Principal  Investigator  was  Dr.  Chang  W.  Sohn.  Dr.  Thomas  J.  Har- 
tranft  is  Chief,  CEERD-CF-E,  and  L.  Michael  Golish  is  Chief,  CEERD-CF. 
The  associated  Technical  Director  was  Dr.  Paul  A.  Howdyshell,  CEERD- 
CF-F.  The  Acting  Director  of  CERL  is  Dr.  Ilker  R.  Adiguzel. 

CERL  is  an  element  of  the  U.S.  Army  Engineer  Research  and  Development 
Center  (ERDC),  U.S.  Army  Corps  of  Engineers.  The  Commander  and  Ex¬ 
ecutive  Director  of  ERDC  is  COL  James  R.  Rowan,  and  the  Director  of 
ERDC  is  Dr.  James  R.  Houston. 
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1  Introduction 

Background 

Different  release  models  can  yield  significantly  different  dynamic  concen¬ 
tration  profiles  in  a  room  depending  on  the  release  rates  chosen.  Results 
from  two  different  cases  demonstrated  significantly  different  concentra¬ 
tion  profiles  in  the  room  of  interest.  There  is  a  need  to  critically  review  the 
current  model  to: 

1.  Identify  limitation  of  the  zero-th  order  model  (zero-th  order  model: 
c  =  M/V  =  uniform  constant) 

2.  Establish  time  and  length  scales,  identify  field  variables  of  significance 
based  on  typical  HVAC  operating  conditions. 

This  ERDC  Special  Report  includes  documents  several  pieces  of  analysis 
accomplished  in  the  first  year  of  (what  was  planned  to  be)  a  3-year  effort. 
At  ERDC  management’s  direction,  the  project  was  stopped  after  comple¬ 
tion  of  the  following  initial  steps: 

1.  Project  Summary  Presentation  (p  3) 

a.  Introduction 

b.  Synopsis  of  FY05  activities 

c.  Benefits  to  (other)  projects 

d.  Ideas  for  follow-on  work 

2.  Experimental  and  numerical  comparisons 

3.  Experimental  tests  (Appendix  A). 

Objective 

The  overall  objectives  of  this  project  were  to: 

1.  Critically  review  the  current  model 

2.  Formulate  a  new,  1st  order  Algebraic  model 

3.  Use  Experimental  data  to  validate  the  modeled  theory. 

Approach 


Work  in  this  project  included: 

1.  Literature  Review 

2.  Initial  Research  Results 
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3.  Collaboration  Development  with  Lawrence  Berkeley  National  Labora¬ 
tory  [LBNL]  Experimental  Data) 

Scope 

This  report  documents  preliminary  work  that  was  suspended  after  Fiscal 
year  2005  (FY05).  No  follow-on  research  is  currently  funded. 
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2  Project  Summary 


ly  Engineer  Research  &  Development  Center  IS  I  am  US  Army  Engineer  Research  &  Development  Center 
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1.  Introduction 


-Problem  Statement 

-FY05  Scheduled  Activities 

-Impact  of  Source  BC 
on  dispersion  modeling 


One  Corps,  One  Regiment,  One  Team  . . .  Serving  Soldiers,  the  Army,  the  Nation 

Problem  Statement 


•  Formulate  concentration  field  at  the  initial  release  of 
particles  from  source  (0th,  Algebraic,  and  differential 
models) 


/ 

*  / 

o  o  o  q.  Particle 

°  1  ®°!00o°Wo  0 concentration 

- o°  ••  .•  o  •  •  •  o  °  ®  O  o 

O  a  O  ri  <Po  O  °  o  o °  °  O  O  O  ° p  O  n 

c(0= 


where: 


Particle  material  property  (P)  -»  pp,  a,  r,  d 
Surface  property  (S) ->  R,  g,... 

Air  stream  property  (A)  ->  v  ,  T,  pa  p,  h 


•  HYPOTHESIS:  Hiaher  order  models  provide  realistic  particle 


One  Corps,  One  Regiment,  One  Team  . .  .  Serving  Soldiers,  the  Army,  the  Nation 


US  Army  Engineer  Research  &  Development  Center  gj)  US  Army  Engineer  Research  &  Development  Center 
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Scheduled  Activities  (FY05) 


1.  Critical  review  of  current  model 

a.  Identify  limitation  of  the  Oth  order  model  (0th 
order  model:  c  =  M/V  =  uniform  constant) 

b.  Establish  time  and  length  scales,  identify  field 
variables  of  significance  based  on  typical  HVAC 
operating  conditions 


2.  Algebraic  formulation  (1st  order  model) 

a.  Identify  app  dimensionless  groups 

b.  Dimensional  analysis,  algebraic  formulation 

c.  Dispersion  simulation  on  PAR3D  (CFD  platform) 


3.  Experimental  data  for  theory  validation 
a.  Laboratory  exp  data  from  LBNL 


f  One  Corps ,  One  Regiment,  One  Team  .  . .  Serving  Soidiers,  the  Army,  the  Nation 


•  Resuspension  models  (e.g.,  c0  =  const) 
result  in  significant  differences  in  the 
dynamic  concentration  profiles  in  a  room 
depending  on  the  values  chosen. 

•  Results  from  two  different  cases  (c0  =  i 

gram/sec  for  10  minutes,  and  C0  =  12  gram/ms  for  50 

milliseconds)  demonstrated  significantly 
different  concentration  profiles  in  the 
room  of  interest. 

f  One  Corps,  One  Regiment,  One  Team _ Serving  Soidiers,  the  Army,  the  Nation 
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600  gram  in  10 
minutes  and  50 
milliseconds: 

After  1  min 


IMPACT  OF  RELEASE  MECHANISM  ON  SIMULATION 


SLOW 


600  gm  released  continuously  at  1  gm/ 

CONCENTRATION  AT  1  MIN 


7.19  mg/I 


600  gm  released  impulsively  at  time  zero 


4  FT  FROM  FLOOR  -  UPPER  GROUND  LEVEL 


IMPACT  OF  RELEASE  MECHANISM  ON  SIMULATION 


600  gram  in  10 
minutes  and  50 
milliseconds: 

After  10  min 


SLOW 


600  gm  released  continuously  at  1  gm/sec 

CONCENTRATION  AT  10  MIN 


FAST 


600  gm  released  impulsively  at  time  zero 


4  FT  FROM  FLOOR  -  UPPER  GROUND  LEVEL 


US  Army  Engineer  Research  &  Development  Center  IfflM  US  Army  Engineer  Research  &  Development  Center 
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Identification  of  scales  of  field  variables 


Definition  ofscafes  ofjnterest 

Length  scale 

-  Room  dimension  in  m  (Lr) 

-  Particle  dimension  in  jn  (dp) 

-  Roughness  dimension  mm  (R) 

-  Boundary  layer  thickness  in  Cm  (8),  based  on  Blasius  solution 
Mass  scale 

-  Particle  density  in  gram/Cm3  (pp) 

-  Air  density  in  0.0012  gram/Cm3  (pa) 

Time  scale  in  second 
Velocity  scale 

-  Room  air  bulk  motion  in  m/sec  (ur) 

-  Particle  settling  velocity  in  mm/sec  (up),  based  on  Stokes  flow 

-  Friction  velocity  in  m/sec  (u„) 

Concentration  scale  in  ppb  (c) 


One  Corps,  One  Regiment,  One  Team  . . .  Serving  Soldiers,  the  Army,  the  Nation 


US  Army  Engineer  Research  &  Development  Center  \  n  US  Army  Engineer  Research  &  Development  Center 
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Accomplishments  -  Literature  Review 


Literature  review  in  the  related  areas: 

-  Indoor  Air  Quality  Modeling 

-  Outdoor  Resuspension  Modeling 

-  Erosion  and  Sedimentation  Modeling 


(Slides  on  each  topic  are  attached  at  the 
end  of  presentation.) 


One  Corps,  One  Regiment,  One  Team  . . .  Serving  Soldiers,  the  Army,  the  Nation 


Accomplishments  -  Algebraic  Model 


Algebraic  model  development: 
C  (t)  =  f  (Re*) 


Preliminary  linear  algebraic  model: 

C  (tV  Cp  =  a  Re*m/Re*„ 


-  Note  that  Re*  (t)  is  a  function  of  time  due  to  u*(t). 

-  Re*c  is  the  critical  Reynolds  number  for 

resuspension 

-  C0  is  specified  by  initial  flux  density,  i.e.,  C0  =  M/S 

-  a  is  an  empirical  constant  determined  from  data. 


One  Corps,  One  Regiment,  One  Team  . . .  Serving  Soldiers,  the  Army,  the  Nation 


US  Army  Engineer  Research  &  Development  Center 
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Accomplishments  -  Experimental  Validation 


Collaboration  with  LBNL 


scoping  experiments  will  be  performed  to  determine  the 
potential  significance  of  resuspension  caused  by  turbulence 
in  the  indoor  environment.  The  basic  experimental 
procedure  will  involve  the  deposition  of  a  known  quantity  of 
a  well  characterized  poly-dispersed  particle  source  (such  as 
ISO  12103-1,  A1  ultra-fine  test  dust)  onto  a  surface, 
subjecting  the  test  surface  to  well  characterized  turbulent 
flow  inside  an  experimental  chamber,  and  measuring  the 
airborne  concentration  of  resuspended  particles  within  the 
chamber  using  a  particle  sizing  instrument  (TSI,  Inc.:  APS 
3321)  to  determine  the  size  dependent  resuspension  rate.“ 


One  Corps,  One  Regiment,  One  Team _ Serving  So/diers,  the  Army,  the  Nation 


V  One  Corps,  One  Regiment,  One  Team _ Serving  Soldiers,  the  Army,  the  Nation 


LBNL  Test  Chamber 


Test  chamber  used  in  the  LBNL  experiments 


Fihiii 


initrtl 


Samples 


US  Army  Engineer  Research  &  Development  Center  IS  J :HK:  US  Army  Engineer  Research  &  Development  Center 
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PRODUCTS  in  FY05 


1.  White  Paper,  “White  Paper  on  Gas/Particle  Dispersion 
in  Armv  Facilities  (AT23-05^2),  “  11  JAN  05 


2.  PAR3D  Validation,  “Flow  and  Dispersion  in  a  Pseudo- 
Rectangular  Chamber.  Comparison  of  PAR3D 
Predictions  with  Experiments  Conducted  by  Shimada 
et.  al.,”  6  Apr  05. 

3.  LBNL  Experimental  Data,  “Experimental 
Measurements  of  Particle  Resuspension  by  Airflow  in 
a  Bench-Scale  Chamber:  Preliminary  Results."  Sep 
05 

4.  Project  Wrap  UP  Briefing  PP  file,  “Resuspension 
Physics  of  Fine  Particles.”  17  Oct  05. 

(These  products  are  provided  in  electronic  files.) 

ft™  J 

One  Corps,  One  Regiment,  One  Team  . . .  Serving  Soldiers,  the  Army,  the  Nation 


3.  Benefits  to  Current  Projects 


Product  #2,  “PAR3D  Validation  -  Flow  and 
Dispersion  in  a  Pseudo-Rectangular 
Chamber.  Comparison  of  PAR3D 
Predictions  with  Experiments  Conducted  by 
Shimada  et.  al.*”  is  used  as  a  part  of  6.2  M&S 
CFD  platform  Verification  and  Validation. 


[*  Shimada,  M.,  K.  Okuyama,  S.  Okazaki,  T.  Asai,  M. 
Matsukura,  and  Y.  Ishizu,  “Numerical  Simulation  and 
Experiment  on  the  Transport  of  Fine  Particles  in  a 
Ventilated  Room,”  Aerosol  Science  and  Technology. 
Vol  25,  Issue  3,  October  1996,  pp.  242-253] 


One  Corps,  One  Regiment,  One  Team  . . .  Serving  Sofdiers,  the  Army,  the  Nation 
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4.  Ideas  for  Follow  On  Works 


1.  Data  Reduction  of  the  LBNL  Experiments 
in  the  form  of  c  (t)  =  f  (Re*) 


2.  Review  of  the  algebraic  model  (C  (t)/  Cn  = 
a  Re*(t)/Re*J  in  comparison  to  the  LBNL 
experimental  data 


3.  Refinement  of  the  algebraic  model  for 
application  to  general  dispersion 
transport  modeling  (multizone,  CFD,  etc.) 


One  Corps,  One  Regiment,  One  Team  . . .  Serving  Soldiers,  the  Army,  the  Nation 


1.  In  ERDC  -  Another  try  in  6.1  competition? 


2.  In  DoD  -  Look  for  basic  research 
opportunities. 


3.  In  Academia  -  Team  up  with  academia 
(e.g.,  UIUC)  for  a  proposal  to  NSF? 


4.  In  Public  -  A  Research  Topic  for  ASHRAE 
TC  4.10  (Indoor  Environmental 
Simulation) 

One  Corps,  One  Regiment,  One  Team _ Serving  Soidiers,  the  Army,  the  Nation 


US  Army  Engineer  Research  &  Development  Center  {  im  nil  US  Army  Engineer  Research  &  Development  Center 
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Additional  Slides  for 
Detailed  Information 


One  Corps,  One  Regiment,  One  Team _ Serving  Soldiers,  the  Army,  the  Nation 
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Two  Kev  Models  in  M&S: 


D 

osaae  an 

a 

source  mo 

5HJ 

Prediction  accuracy  of  concentration  can  be  only  as  accurate 
as  the  source  strength  input. 


No  accepted  theoretical  model  of  indoor  resuspension  vet : 
(Subject  of  this  Project) 


Experiments  currently  on  going  for  empirical  model  development 
(LBNL,  LLNL) 

Dosage  =  Concentration  *  Exposure  time*  Retention  efficiency 
(Time-varying  concentration  at  a  point  is  the  Life  Safety  Issue.) 
-Beyond  the  scope  of  this  project. 


“Source  term  models  and  dose  response  models  are  the  two  key 

technical  elements”  for  “existing  and  future  interior  and 
exterior  codes.”  [SNL/ECBC,  DHS  Conf  paper,  Apr  05] 


f  One  Corps ,  One  Regiment,  One  Team  .  . .  Serving  Soidiers,  the  Army,  the  Nation 


•  Limitation  of  the  0th  order  model  in 
dynamic  simulation  of  concentration 
distribution  in  a  room  is  demonstrated  - 
The  inherent  assumption  of  “uniform 
concentration  in  a  space,  i.e. ,  c=M/V” 
cannot  provide  desired  concentration 
information  in  a  room  of  interest. 


One  Corps,  One  Regiment,  One  Team  . .  .  Serving  Soldiers,  the  Army,  the  Nation 
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Accomplishment  -  Limitation  of  Oth  order  model 


r  One  Corps,  One  Regiment,  One  Team _ Serving  Soidiers,  the  Army,  the  Nation 


Accomplishment  -  Algebraic  Model 


Identified  Dimensionless  groups 

•  Red  =  upd/i>  (Particle  Reynolds  number) 

•  Res  =  ur  5  !v  (Re  based  on  turbulent  boundary  layer 
thickness,  8) 

•  Fr  =  u*2/gS  (Froude  number) 

•  dp/R  =  (Roughness  factor) 

•  St  =  rp/rL(Stokes  number) 

•  Gr  =  p2gd2Ap/ji2  (Grashof  number) 

•  DF  =  upoo/u*  (drift  factor) 

[“Numerical  approaches  for  motion  of  dispersed  particles,  droplets  and  bubbles,”  E. 
Loth,  2000,  Progress  in  Energy  and  Combustion  Science,  vol.  26,  pp.  161-223] 


One  Corps,  One  Regiment,  One  Team _ Serving  Soldiers,  the  Army,  the  Nation 
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Accomplishments  -  Literature  Review 


Related  Studies  on  Indoor  Air  Quality  (IAQ)  Modelin 


In  IAQ  modeling,  either  by  multizone  or  CFD  models,  ...  - 
’’Contaminants  introduced  at  the  inlets  can  be  specified 
as  either  constant  or  varying  according  to  user-defined 

functions”  “An  Analysis  of  Combined  CFD  and  Multizone  IAQ  Model  Assembly  Issues,” 
A.  Musser,  ASHRAE  Trans,  AT-01-1-3,  pp.  371-382,  2001 


“Choosing  a  Multiphase  Model:  User  must  know  a  priori  the 
characteristics  of  flow:  Flow  regime,  Laminar/turbulent, 
Dilute/dense,  Particulate  Loadinq.Stokes  number, 
Secondary  phase  (??)  diameter  for  drag  considerations” 


“FLUENT  v6.Q:  Training  Notes,  pp.  9-8,  Jan  2002” 


One  Corps ,  One  Regiment,  One  Team  . . .  Serving  Soldiers,  the  Army,  the  Nation 


Accomplishments  -  Literature  Review 


Related  Studies  on  Outdoor  Resuspension 

Resuspension  rate  (A)  and  Resuspension  factor  (K): 

A  =  fraction  of  surface  deposits  removed  per  unit  time 

K  =  (airborne  concentration  of  deposited  material)  /  (quantity  of 
deposited  material  per  unit  area  (projected)  ground  surface. 

Experimental  determination  of  A  from  wind  tunnel  data 
does  not  provide  onset  criteria  for  particle  resuspension. 

[Ref:  “Evaluation  and  development  of  models  for  resuspension  of  aerosols  at  short  times 
after  deposition,”  G.A.  Loosemore,  Atmospheric  Environment,  vol  37,  2003,  pp. 
639-647] 

[Ref:  “Wind  tunnel  experiments  on  the  resuspesnion  of  particulate  material,”  K.W. 
Nicholson,  Atmospheric  Environment,  vol  27A,  1993,  pp.  181-188] 


f  One  Corps,  One  Regiment,  One  Team _ Serving  Soldiers,  the  Army,  the  Nation 


US  Army  Engineer  Research  &  Development  Center  l  ilc  US  Army  Engineer  Research  &  Development  Center 
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Accomplishments  -  Literature  Review 

Outdoor  Empirical  Models  (Current  state  of  art): 

A  =  [u„c1dpc2]  /  [tc3z0c4rpc5]  ,  where 

A  :  Resuspension  rate  (Resuspension  flux/initial  concentration) 

u*:friction  velocity,  dp:  particle  diameter,  t :  time,  z0:  surface 

roughness,  rp:  particle  density,  cl . c5:  empirical 

constants 

Current  models  are  based  on  field  and  wind  tunnel  experimental 
data  for  soil,  grass  and  concrete  surfaces  under  typical 
atmospheric  environment  (Loosmore,  2003). 

e.g.,  A  =  0.42*(u,2 13dp°  17/t?  92z0°  32pp°  76)  [s_1]  (Loosmore,  2003) 


^  One  Corps,  One  Regiment,  One  Team _ Serving  Soldiers,  the  Army,  the  Nation 


'%  Accomplishments  -  Literature  Review 

Related  Studies  on  Erosion  and  Sedimentation 

Resuspension  Condition:  Beginning  of  motion  depends  on  whether  the  boundary 
layer  shear  stress  (xq)  becomes  equal  to  or  greater  than  the  critical 
boundary  layer  shear  stress  (xc) . 

x*c  =  V  xc .  =  f  (Re*);  Re*  =  u*  ds/vm 

Shields  parameter  (dimensionless  shear  stress,  x*  ): 

**  =  V(Ys^m)ds  =  Pmu*2  /(Ys^m)ds,;  where, 

XqI  boundary  layer  shear  stress 
u*:  shear  velocity 

ys:  specific  weight  of  a  sediment  particle 
ym:  specific  weight  of  the  fluid  mixture 
ds :  particle  size 

[Ref:  "Fluid-Particle  Interactions  and  Resuspension  in  Simple  Shear  Flow,"  Z.  Feng  and  E.  Michaelides,  J.  Hydraulic  Eng,  ASCE,  Dec  2003, 
pp.  985-994] 

[Ref:  Erosion  and  Sedimentation.  P.  Julien,  Cambridge  University  Press,  1995,  pp.  114] 

Grain  Shear  Reynolds  number  Re*  (=  u*  ds/vm)  is  the  deciding  factor  in  the 
initiation  of  the  particle  motion. 

One  Corps,  One  Regiment,  One  Team  . . .  Serving  Soldiers,  the  Army,  the  Nation 
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LBNL  Results  -  flow  field  @Fans  at  110 


Velocities  at  each  measuring  point 


7 

8 

9 

Min 

0.35 

0.15 

0.14 

Ave 

0.48 

0.19 

0.27 

Max 

0.74 

0.24 

0.47 

4 

5 

6 

Min 

0.40 

0.22 

0.16 

Ave 

0.57 

0.32 

0.25 

Max 

0.72 

0.41 

0.34 

i 

2 

3 

Min 

0.34 

0.26 

0.20 

Ave 

0.49 

0.49 

0.32 

Max 

0.82 

0.67 

0.46 
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(Provisional)  Project  Chronolo 


-  1  Oct  04:  Project  Start 

-19  Oct  04:  TD  e-mail  on  IPR  prep  “..focus  on  the  science  accomplished  to 
obtain  new  knowledge  and  not  on  application  of  the  knowledge.” 

-  28  Oct  04,  BREP  IPR:  (Experimental  validation  required.) 

-  14  Jan  05:  White  paper  produced  per  comments  from  BREP  IPR 

-  28  Feb  05:  2-page  Quad  Charts  briefing  to  ERDC  PMB 

-  21  Mar  05:  PMB  de-brief 

-  12  May  05:  Re-briefing  to  PMB  (Relevant  to  the  Army?) 

-  12  Jul  05:  6.1  VTC  scheduled/cancelled.  This  project  will  end  in  FY05. 

Wrap  up  FY05  activities  in  a  tech  report. 

-  3  Aug  05:  E-mail  notice,  FY06  Basic  Research  IPR  on  13-16  Sep  05 
- 19  Aug  05:  E-mail  (from  Dr.  Cullinane),  Project  Wrap  UP  VTC 
-17  Oct  05:  Wrap  Up  briefing  to  BREP  TD 


One  Corps,  One  Regiment,  One  Team  . . .  Serving  So/diers,  the  Army,  the  Nation 

(Provisional)  On  PMB  Comments  in  May  05 


“Why  weaponized  anthrax  was  de-ionized?” 


William  C.  Patrick  III,  a  microbiologist*  who  designed  germ  weapons  for  the  United 
States  before  President  Richard  M.  Nixon  renounced  them  in  1969,  said  he  had 
learned  details  of  the  federal  inquiry  from  a  senior  investigator.  The  Senate  powder, 
Mr.  Patrick  said,  was  quite  potent  and  capable  of  sailing  far  through  the  air  to  hurt 
many  people. 

He  said  the  makers  of  the  anthrax  spores  sent  to  Mr.  Daschle's  office  had  produced  a 
dry  powder  that  was  remarkably  free  of  extraneous  material. 

"It  s  high-grade,"  said  Mr.  Patrick,  who  consults  widely  on  making  germ  defenses.  "It's 
free  flowing.  It's  electrostatic  free.  And  it's  in  high  concentration. 

Experts  on  germ  weaponry  agree  that  the  removal  of  electrostatic  charges  is  a  major 
step  toward  making  an  effective  munition.  The  Soviet  Union  and  United  States 
developed  sophisticated  ways  of  diminishing  this  attraction  and  helping  the  particles 
float  more  freely,  increasing  their  ease  of  dissemination  and  likelihood  of  inhalation. 

Mr.  Patrick  said  that  whoever  sent  the  Daschle  letter  had  clearly  achieved  this  step. 
"It's  fluffy  "  he  said,  quoting  experts  who  examined  the  powder.  "It  appears  to  have  an 
additive  that  keeps  the  spores  from  clumping."  Removing  the  charge,  he  added,  is  a 
black  art,  few  details  of  which  are  known  publicly. 

Assertions  by  some  federal  officials  that  the  material  was  not  the  type  that  would  be 
used  in  weapons  are  "nonsense,"  he  said.  "The  only  difference  between  this  and 
weapons  grade  is  the  size  of  the  production.  ...” 

[The  New  York  Times,  October  25,  2001 ,  “Contradicting  Some  U.S.  Officials,  3 
Scientists  Call  Anthrax  Poweder  High-Grade] 

U  *William  C.  Patrick  III,  a  microbiologist  -  Former  Director  of  Operation  at  Ft.  Detrick 

y  One  Corps,  One  Regiment,  One  Team _ Serving  So/diers,  the  Army,  the  Nation 
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(Provisional)  On  PMB  Comments  in  May  05 


“Place  fine  material  simulating  anthrax  on  an 
accurate  scale.  Leave  it  alone  over  the  weekend  and 
see  how  much  remains  when  you  come  back.  Even 
such  a  back-of-the  envelope  approach  has  not  been 
done  to  see  if  the  effort  would  have  even  slight 
utility .” 

LBNL  experimental  data  in  the  Product  #3  provide  the 
answers  to  the  comment. 
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3  Flow  and  Dispersion  in  a  Pseudo 
Rectangular  Chamber 


Flow  and  Dispersion  in  a 
Pseudo-Rectangular  Chamber 

Comparison  of  PAR3D  Predictions  with 
Experiments  Conducted  by  Shimada  et  al. 

Performed  by  Robert  S.  Bernard,  ERDC-CHL 


Grid  outline  indicates  shape  of  test  chamber 
and  locations  of  inflow  and  outflow  ports. 

Chamber  is  270  cm  wide  (x-direction), 

330  cm  long  (y-direction),  and 
200  cm  high  (z-direction). 

Air  flows  into  chamber  through  slot  (5  cm  wide) 
from  elliptical  supply  duct  along  ceiling. 

Air  flows  out  of  chamber  through 
40-cm  square  opening  centered  at 
x  =  135  cm,  y  =  0,  and  z  =  35  cm. 

Airflow  rate  =  257  liters  per  second 

Only  steady-state  results  are  shown  for  velocity. 
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Velocity  Vectors  in  Central  YZ- Plane 

Computed  with  no-slip  condition 


Velocity  Vectors  in  Central  YZ- Plane 
Computed  with  slip  condition 
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Velocity  Vectors  in  Central  XZ- Plane 

Computed  with  no-slip  condition 
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Velocity  Vectors  in  Central  XZ- Plane 

Computed  witli  slip  condition 


z-component  Of  velocity.  m/S  Z-component  of  velocity,  m/s 
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Experimental  data  for  z-component  of  velocity  in  central  xz-plane  (•)  compared 
with  PAR3D  predictions  using  slip  condition  ( - )  and  no-slip  condition  ( - ) 


0  0.3  0.6  0.9  1.2  1.5  1.8  2.1  2.4  2. 


Y-position.  m 


Experimental  data  for  z-component  of  velocity  in  central  xz-plane  (•)  compared 
with  PAR3D  predictions  using  slip  condition  ( - )  and  no-slip  condition  ( - ) 
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0  0.3  0.6  0.9  1.2  1.5  1.8  2.1  2.4  2.7 

Y-position.  m 


Experimental  data  for  z-component  of  velocity  in  central  xz-plane  (•)  compared 
with  PAR3D  predictions  using  slip  condition  ( - )  and  no-slip  condition  ( - ) 


A  particle  emitter  is  placed  at  the  intersection 
of  the  central  xz-  and  yz-planes, 
either  50  cm  above  the  floor  (emitter  A) 
or  directly  on  the  floor  (emitter  B). 

The  normalized  emission  rate  is 
0.0167  units  per  second, 
equivalent  to  a  concentration  of 
1.0  units  per  liter  at  an  airflow  rate 
of  0.0167  liters  per  second. 

Particle  concentrations  shown  here 
are  concentrations  relative 
to  the  concentration  at  the  emitter. 

Only  steady-state  results  are  shown 
for  particle  concentration. 


ERDC/CERL  SR-06-18 


29 


Relative  Concentration  in  Central  YZ-Plane  for  Emitter  A 

Computed  with  no-slip  condition 
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Relative  Concentration  in  Central  YZ-Plane  for  Emitter  A 

Computed  with  slip  condition 


Relative  Concentration  in  Central  YZ-Plane  for  Emitter  B 

Computed  with  no-slip  condition 
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Relative  Concentration  in  Central  YZ-Plane  for  Emitter  B 

Computed  with  slip  condition 


Relative  Concentration  in  Central  XZ- Plane  for  Emitter  A 

Computed  with  no-slip  condition 


ERDC/CERL  SR-06-18 


32 


ERDC/CERL  SR-06-18 


33 


Relative  Concentration  in  Central  XZ-Plane  for  Emitter  B 

Computed  with  slip  condition 


Experimental  data  for  relative  particle  concentration  in  central  xz-plane  (•)  compared 
with  PAR3D  predictions  using  slip  condition  ( - )  and  no-slip  condition  ( - ) 
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4  Summary 


A  chronological  summary  of  this  research  activity  follows: 


1.  (l  Oct  04) 

2.  (19  Oct  04) 


3.  (28  Oct  04) 

4.  (14  Jan  05) 

5.  (28  Feb  05) 

6.  (21  Mar  05) 

7.  (12  May  05) 

8.  (12JUI05) 

9.  (3  Aug  05) 

10.  (19  Aug  05) 

11.  (17  Oct  05) 


Project  Start 

TD  e-mail  on  IPR  prep  “..focus  on  the  science  accom¬ 
plished  to  obtain  new  knowledge  and  not  on  applica¬ 
tion  of  the  knowledge.” 

BREP  IPR:  (Experimental  validation  required.) 

White  paper  produced  per  comments  from  BREP  IPR 
2-page  Quad  Charts  briefing  to  ERDC  PMB 
PMB  de-brief 

Re-briefing  to  PMB  (Relevant  to  the  Army?) 

ERDC  Direction:  This  project  will  end  in  FY05. 

(Wrap  up  FY05  activities  in  formal  report.) 

E-mail  notice,  FY06  Basic  Research  IPR  on  13-16  Sep 

05 

E-mail  (from  Dr.  Cullinane),  Project  Wrap  UP  VTC 
Wrap  Up  briefing  to  BREP  TD. 
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Appendix  A:  Experimental  Measurement  of 
Particle  Re-Suspension  by  Airflow  in  a 
Bench-Scale  Chamber:  Preliminary 
Results 

Mark  Sippola 

Indoor  Environment  Department,  Lawrence  Berkeley  National  Laboratory 

Introduction 

Particles  that  settle  onto  floors  of  indoor  environments  can  re-suspend 
into  the  air  and  be  transported  to  other  locations.  Particles  are  known  to 
be  re-suspended  by  human  activities,  but  the  extent  to  which  particles  may 
be  re-suspended  by  airflows  indoors  is  unclear.  A  set  of  exploratory  ex¬ 
periments  have  been  conducted  to  better  understand  the  role  that  airflows 
may  play  in  re-suspending  particles  settles  on  floor  surfaces  indoors. 

Methods 

Particle  re-suspension  by  airflows  was  measured  in  a  series  of  four  ex¬ 
periments.  The  experiments  investigated  the  re-suspension  of  particles  in 
the  size  range  1.7-3. 8  and  5.0-10.3  pm  from  floors  covered  with  lino¬ 

leum  and  carpet.  The  conditions  of  the  experiments  are  summarized  in 
Table  1. 

The  experiments  were  conducted  in  a  30x30x36  inch  (LxWxH)  chamber 
with  two  variable  speed  fans  mounted  on  the  interior  to  induce  airflow 
within  the  chamber.  The  floor  of  the  chamber  was  completely  covered  with 
either  moderately  textured  linoleum  or  commercial  grade  olefin  carpet. 
Please  refer  to  the  MS  Word  file  Re-suspension  -  Physical  description  for 
a  more  detailed  description  of  the  chamber.  Each  experiment  had  four 
phases:  pre-cleaning,  deposition  of  particles  onto  floor  surfaces,  re¬ 
suspension  of  particles,  and  quantification  of  particle  mass  on  floor  sam¬ 
ples.  Please  refer  to  the  MS  Word  file  Re-suspension  -  General  procedure 
for  a  more  detailed  description  of  the  experimental  procedure. 

During  the  particle  deposition  phase,  the  airborne  particle  concentration 
was  measured  with  an  Aerodynamic  Particle  Sizer  (APS)  (TSI,  Inc.,  Model 
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3321).  The  APS  measured  airborne  particle  concentration  in  51  size  bins  in 
the  range  0.5-20  |Lim.  Data  collected  by  the  APS  were  used  to  estimate  the 
size  distribution  of  particles  deposited  on  the  floor  surfaces. 

During  the  particle  re-suspension  phase,  the  airborne  concentration  was 
measured  with  both  an  APS  and  a  filter  sample.  The  re-suspension  phase 
for  each  experiment  lasted  5  hours,  1  hour  with  the  fans  off  and  1  hour  for 
each  of  four  fan  speeds  (45,  75,  no  and  140).  The  fan  speeds  were  incre¬ 
mentally  increased  each  hour  over  the  course  of  the  experiment.  Filters 
samples  were  collected  for  each  hour  at  a  given  fan  speed.  The  APS  con¬ 
tinuously  collected  data  which  was  integrated  over  10  minute  time  inter¬ 
vals. 

Air  velocities  in  the  chamber  were  measured  with  a  hotwire  anemometer 
(Airflow  TA5)  for  all  5  fan  speeds.  Measurements  were  made  at  9  points  in 
a  plane  10  cm  above  the  floor.  The  9  measurement  points  were  at  the  cen¬ 
ter  points  of  a  9-square  grid  in  the  plane  parallel  to  the  floor.  Measure¬ 
ments  were  made  with  the  probe  oriented  in  both  the  horizontal  and  verti¬ 
cal  directions.  The  hotwire  anemometer  recorded  a  measurement  about  3 
times  per  second.  It  was  held  at  each  measurement  location  for  30  seconds 
and  then  the  minimum,  average  and  maximum  air  velocity  in  that  30  sec¬ 
ond  time  frame  was  recorded. 

There  are  two  potential  weaknesses  of  these  experiments.  During  the 
deposition  phase,  some  particles  may  deposit  to  the  walls  and  other  sur¬ 
faces  within  the  chamber.  These  particles  that  are  not  deposited  on  the 
floor  surface  are  potentially  available  for  re-suspension.  The  fans  lead  to 
some  vibration  of  the  chamber  when  they  were  turned  on,  and  the  inten¬ 
sity  of  the  vibration  increased  as  the  fan  speed  was  increased.  This  vibra¬ 
tion  could  potentially  lead  to  particle  re-suspension  independent  of  the 
airflows.  Further  investigation  will  be  required  to  understand  the  degree  to 
which  two  these  factors  may  have  influenced  the  experimental  results. 

Results 

The  air  velocities  measured  in  the  chamber  for  the  5  different  fan  speeds 
are  presented  in  Figure  1.  The  values  presented  are  the  average  values  for 
the  nine  air  velocity  measurement  locations.  The  inflow  of  make-up  air 
and  the  outflow  of  sampled  air  to  the  APS  and  filter  lead  to  the  nonzero  air 
velocities  measured  when  the  fans  are  off.  The  air  velocities  at  fan  speeds 
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of  75,  no  and  140  are  similar.  The  fan  speeds  were  originally  selected 
based  on  the  volumetric  flow  rates  presented  in  Table  2. 

Particle  re-suspension  was  not  observed  in  any  experiment  when  the  fans 
were  off  or  at  speed  45.  In  all  four  experiments,  particle  re-suspension  was 
observed  at  fan  speeds  of  75, 110  and  140.  Figure  2  shows  a  series  of  plots 
of  airborne  concentration  versus  particle  size  data  collected  by  the  APS 
during  the  re-suspension  phase  of  Experiment  2.  Each  plot  in  the  figure 
shows  airborne  particle  concentration  integrated  over  successive 
10-minute  intervals  beginning  at  the  point  when  the  fans  were  increased 
from  speed  45  to  75.  This  series  of  data  displays  trends  that  were  also  ob¬ 
served  at  higher  fan  speeds  and  during  different  experiments.  Within  a  few 
seconds  after  increasing  the  fan  speed  to  75,  a  large  increase  in  the  air¬ 
borne  particle  concentration  was  detected  by  the  APS  in  the  size  range  1-5 
pm.  This  concentration  remained  fairly  constant  for  about  10  minutes  and 
then  began  to  decrease.  Figures  2a  through  2e  show  a  decrease  in  airborne 
particle  concentrations  while  the  fan  speed  remains  at  75. 

Figure  3a  shows  the  fraction  of  deposited  particle  mass  that  was  re¬ 
suspended  from  the  linoleum  floor  in  1  hour  while  the  fans  were  at  a  speed 
of  75  for  a  range  of  particle  sizes  and  Figure  3b  shows  the  same  informa¬ 
tion  for  the  carpeted  floor.  The  data  in  Figure  3a  come  from  combining  the 
data  collected  for  smaller  particles  in  Experiment  2  with  data  collected  for 
larger  particles  in  Experiment  4.  The  discontinuity  in  the  data  at  a  particle 
size  of  about  4  pm  is  where  the  two  data  sets  meet  but  disagree.  The  data 
shown  in  Figure  3b  come  from  small  particle  data  collected  in  Experiment 
1  and  large  particle  data  collected  in  Experiment  3. 

Figures  4  and  5  show  the  same  data  as  Figure  3  for  fan  speeds  of  no  and 
140,  respectively.  In  these  figures,  the  data  for  small  and  large  particles 
agree  and  there  is  no  strong  discontinuity  at  the  transition  between  the 
two  data  sets. 


Table  Ai.  Summary  of  experimental  conditions.. 


Experiment  #  (-) 

Floor  surface  (-) 

Particle  size  range  (|jm) 

Initial  mass  lading  (ng/cm2) 

1 

carpet 

1-4 

460 

2 

linoleum 

1-4 

160 

3 

carpet 

5-10 

5200 

4 

linoleum 

5-10 

6700 
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Figure  Ai.  The  average  of  the  minimum,  average  and  maximum  air  velocity 
in  the  (a)  horizontal  and  (b)  vertical  direction  measured  at  the  nine 
measurement  locations  for  the  five  different  fan  speeds. 
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(e) 

Figure  A2.  Airborne  mass  concentration  of  particles  measured  by  the  APS 
during  the  re-suspension  phase  of  Experiment  2.  The  data  are  10-minute 
integrated  concentrations  collected  during  the  following  time  intervals  after 
the  fan  speed  was  increased  to  75:  (a)  0-10  minutes,  (b)  10-20  minutes, 

(c)  20-30  minutes,  (d)  30-40  minutes  and  (e)  40-50  minutes. 
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Figure  A3.  Fraction  of  deposited  particle  mass  re-suspended  from  (a) 
linoleum  and  (b)  carpet  during  1  hour  at  a  fan  speed  of  75. 


Figure  A4.  Fraction  of  deposited  particle  mass  re-suspended  from 
(a)  linoleum  and  (b)  carpet  during  1  hour  at  a  fan  speed  of  110. 
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Figure  A5.  Fraction  of  deposited  particle  mass  re-suspended  from 
(a)  linoleum  and  (b)  carpet  during  1  hour  at  a  fan  speed  of  140. 
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